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We consider non-standard cosmological models in which the late decay of a scalar field <j> reheats the 
Universe to a low reheating temperature, between 5 MeV and the standard freeze-out temperature 
of neutralinos of mass m x . We point out that in these models all neutralinos with standard density 
fistd ^ 10~ J (100GeV/m x ) can have the density of cold dark matter, provided the right combination 
of the following two parameters can be achieved in the high energy theory: the reheating temper- 
ature, and the ratio of the number of neutralinos produced per (ft decay over the eft field mass. We 
present the ranges of these parameters where a combination of thermal and non-thermal neutralino 
production leads to the desired density, as functions of fi st d and m x . 

PACS numbers: 14.60.St, 98.80.Cq 



In supersymmetric models, the lightest supersymmet- 
ric particle (LSP), usually a neutralino x, is a good cold 
dark matter candidate. The cosmological density of the 
neutralino is a function of the supersymmetric model 
parameters, and it has been computed theoretically to 
high precision. Requiring the LSP to have the measured 
dark matter density constrains the models considerably 
to very narrow regions in parameter space. 

The standard computation of the relic density relies 
on the assumptions that the entropy of matter and ra- 
diation is conserved, that neutralinos are produced ther- 
mally and were in thermal and chemical equilibrium be- 
fore they decoupled. The decoupling, or freeze-out, tem- 
perature, i.e. the temperature after which their number 
practically does not change, is Tf . ~ m x /20, where 
m x is the neutralino mass. The standard neutralino 
relic density fi st d obtained in this way can be larger or 
smaller than the measured density of cold dark matter 
^ cdm = 0.113 ±0.009/i" 2 0. 

However, there are non-standard cosmological models 
in which the assumptions mentioned above do not hold. 
These include models with moduli decay , Q-ball decay 
, and thermal inflation Q . In all of these models there 
is a late episode of entropy production and non-thermal 
production of the LSP in particle decays is possible. 

We concentrate on cosmological models in which the 
late decay of a scalar field (ft reheats the Universe to a 
low reheating temperature Trh, which can be smaller 
than Tf . without spoiling primordial nucleosynthesis Q 
(5 MeV < Trh ^ Tf.o.)- Such scalar fields are com- 
mon in superstring models where they appear as moduli 
fields. The decay of eft into radiation increases the en- 
tropy, diluting the neutralino number density. The de- 
cay of eft into supersymmetric particles, which eventually 
decay into neutralinos, increases the neutralino number 
density. We denote by b the net number of neutralinos 
produced on average per eft decay. 

Both thermal and non thermal production mechanisms 



have been discussed in the literature 0, IE 0, H, U E3, 
ITTl ll2L Il3| . Several supersymmetric models with par- 
ticular combinations of Tr h b. and m§ have been stud- 
ied 0, H H H [H IH El El Ref. assumed that 
one neutralino is produced per eft decay (b ~ 1), and 
concluded that in the MSSM the non-thermal neutralino 
production leads to an overabundance of relic neutrali- 
nos by several orders of magnitude. Also the thermaliza- 
tion ofy's produced in decays was discussed in Ref. 0. 
Refs. pi lol Il2fl studied the thermal production of LSPs 
during the decay of a scalar field with 6 = 0, and 
mentioned the possibility of b ^= but did not study 
it. Ref. considered purely non-thermal production. 
Ref. 01 is closest to our work, in that both thermal and 
non-thermal production were considered, but the general 
strategy to rescue models with too low or too high stan- 
dard relic density remained, in our opinion, unclear. 

In spite of all the above work, no coherent overview of 
the issue, which is remarkably easy to understand, has to 
our knowledge been laid down. This is what we intend 
to provide in this letter. 

Let J7 s td be the density that neutralinos in a particular 
model would have with the usually assumed cosmology 
An appropriate combination of the following two param- 
eters can bring the relic neutralino density to the desired 
value f2 c dm: the ratio b/m^ giving the average number of 
LSPs produced per unit of energy released in each (ft de- 
cay and the reheating temperature Trh (which must be 
< Tf . ). We assume that the oscillating (ft field dominates 
the energy density of the Universe at early times and that 
at Trh the Universe becomes dominated by the radiation 
produced in (ft decays. We follow the usual choice of the 
parameter Trh as the temperature the Universe would 
attain under the assumption that the (ft decay and sub- 
sequent thermalization are instantaneous 
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F^ is the decay width of the <f> field, 

^ * t/- (2) 

iV eff 

If (f) has non-suppressed gravitational couplings, as is usu- 
ally the case for moduli fields, the effective energy scale 
A c ff ~ Mp. In models with intermediate scales, A c g 
could be smaller [l]J. Thus, with g± ~ 10 

T RH ^10Mev( m * Y /2 (^). (3) 
HH VlOO TcVJ \A cS J y ' 

In order not to disrupt the predictions of Big Bang nu- 
cleosynthesis, Tr H > 5 MeV 0. Thus 

^>100Tev(|g) 2/3 . (4) 

The number b of neutralinos produced per 0-decay is 
highly model-dependent. It is determined by the physics 
of the hidden sector, by the mechanism of supersymme- 
try breaking, and in superstring-inspired models by the 
compactification mechanism. 

The coupling of the <fi to the gravitino arises from the 
term e if ' 2 ^>„<7 A " / ^>„, where K is the Kahler potential. If 
rricj, is larger than twice the gravitino mass m.3/2, the de- 
cay mode <p — > "03/2^3/2 °f the moduli field into two 
gravitinos is present with branching ratio of order 0.01 
(see Refs. Q, which correct previous claims that 
this branching would be chirally suppressed by a factor 
(%/2/™<^) 2 )- Gravitinos must then decay rapidly not 
to disrupt nucleosynthesis (so m 3 / 2 > 100 TeV), and 
they produce comparable amounts of normal particles 
and their supersymmetric partners. If to^ 3> the 
gravitino decays during the radiation dominated epoch 
after the decay of the <f> field (here we do not address 
this case and we focus on neutralino production during p 
domination). When and m.3/2 are of the same order 
of magnitude, we can consider the gravitino decay as part 
of the 4> decay, since they happen almost simultaneously. 
In this case, depending on how important the direct de- 
cay of <f> into supersymmetric particles other than the 
-03/2 is, b can typically be 0.01-1, but not smaller. 

If instead to^ < 2m 3 / 2 more possibilities open up. The 
yield per tfi decay b can still be of order one but it can 
also be much smaller. 

Supergravity models with chiral superfields <E>/ are 
specified in terms of the Kahler potential K(Q], $/), 
the superpotential W($>i), and the gauge kinetic func- 
tion f a p(^i). Specific relations between the <p mass 
m^, the gravitino mass m 3 / 2 , and the gaugino mass 
TO1/2 arise as a consequence of the relations TO3/2 = 
(e K / 2 W), mi/a = (F J djlnRe f), and m<t> = (d 2 V/dcj) 2 ). 
With appropriate choices of K, W, and /, the hier- 
archy TO3/2 > to^ S> mi/2 may be achieved. Here 

V = KjjF 1 ^ -3e K \W\ 2 +^(Ref)-pD a Df ) is the scalar 
potential, F 1 = -e K / 2 K l7 D-jW is the F-term of the 
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FIG. 1: Neutralino density Q, x h 2 as a function of the re- 
heating temperature Trh for illustrative values of the ratio 
77 = ^(lOOTeV/m^) (the number of neutralinos produced per 
unit of energy released in (f> decay, in units of (100 TeV) -1 ). 
The curves are insensitive to the initial conditions, i.e. to the 
value Hi of the Hubble parameter at the beginning of the <j> 
dominated epoch. 

chiral superfield $/, D a = T aIJ ^DjW/W is the D- 
tcrm of the vector superfield, DjF = diT ' + TdjK is the 
Kahler covariant derivative, d\T — dF/d&i, and square 
brackets denote vacuum values. As far as b is concerned, 
one finds, for example, b ~ O(l) when the main <j> de- 
cay mode is through a coupling of the type hipip 2 with 
a chiral matter supermultiplet ip in the superpotential 
W. This leads to comparable decay rates of <fi into the 
scalar and fermionic components of ip (which are super- 
symmetric partners). On the other hand, it is possible 
that the 4> field decays mostly into Higgs fields, or gauge 
fields (Ws, Z's, photons, gluinos). In this case b can 
be very small 10~ 2 , 10~ 4 , 10~ 6 etc. HG3. For exam- 
ple, the coupling of <fi to the gauge bosons arises from 
the term Re f a pF"F ,iV ^ and with non-minimal kinetic 
terms f a p may contain p. The <p decay width into gauge 
bosons is then T g ~ A ff m|/Mj> with A g = ^lnRe/, 
while that into gauginos is Tg ~ X^m^/Mp with Ag = 

m i/2-§^ m (-^^ R- e /)• Thus in principle the gaugino 
coupling may be suppressed relative to the coupling to 
gauge bosons. 

Neutralinos are in kinetic equilibrium with the radi- 
ation when their scattering rate off relativistic particles 
is faster than the Hubble expansion rate, r scatt > H . 
During the epoch in which the Universe is dominated 
by the decaying p> field, H is proportional to T 4 [Tg. 
In fact, in the evolution equation for the radiation en- 
ergy density (Eq. (5) below with p = p/3) substitute 
p ~ T 4 and ~ MpH 2 . Then use H ~ write 
T oc t a , match the powers of t in all terms, and determine 
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a = —(1/4). Hence, during the oscillating cf> dominated 
epoch H oc t^ 1 cx T A (and p^ cx H 2 oc T 8 ). Since F 
equals T| H /M P at T = T RH , it is H ~ T 4 /(T^ H M P ). 
Since r scatt ~ n 7 CT scatt — T 3 a scat t, we see that kinetic 
equilibrium is maintained at low temperatures, T < 
T^jjAfpUgcatt- This is contrary to the usual radiation- 
dominated scenario in which the temperature has to be 
large enough for kinetic equilibrium to be maintained. 
The origin of this difference is the strong T 4 dependence 
of H on the radiation temperature. Using the known 
relation Q, st( ±h 2 — 10 -10 GeV -2 / (av) , and taking <r S catt 
of the same order of magnitude as the annihilation cross 
section (av), gives T < (10 6 MeV /^ st d^ 2 )(T RH /MeV) 2 . 
For the following we need to assume kinetic equilibrium 
before neutralino production ceases: the right hand side 
is larger than T RH for n std h 2 < 10 6 (T RH /MeV), so we 
will safely assume that kinetic equilibrium is reached. 

In kinetic equilibrium the equations which describe the 
evolution of the Universe are: 



p = -3H(p + p) + T <t> p^ 



H 2 = 



-3Hn x - (av) (n\ - n 2 . oq , 



SHpcf, — r^ptp 

— 2 (P + P4>)- 



3M 



(5) 

^<pP4> (6) 
(7) 
(8) 



Notice that the equations, and thus the results, depend 
only on the combination bjm§ and not on b and sep- 
arately. In Eqs. a dot indicates a time deriva- 
tive, pej, is the energy density in the <p field, which is 
assumed to behave like non-relativistic matter; p and p 
are the total energy density and pressure of matter and 
radiation at temperature T, which are assumed to be 
in kinetic but not necessarily chemical equilibrium; n x 
is the number density of LSPs, and ri Xi0q is its value 
in chemical equilibrium; finally, H = a/ a is the Hub- 
ble parameter, with a the scale factor. For convenience 
in the numerical calculations, we used as independent 
variable In a and as dependent variables Y = n x /s, 
Yc/, = p^Km^s) and T, where s = (p + p — m x n x )/T is 
the entropy density of the matter and radiation. We also 
used the first principle of thermodynamics in the form 
d(pa 3 ) + d(p^a 3 ) +pda 3 = Td(sa 3 ) to rewrite Eq. JSJ) as 



s = -SHs 



F<j>p<p 
T 



(9) 



Initial conditions are specified through the value Hj of 
the Hubble parameter at the beginning of the <j) domi- 
nated epoch. This amounts to giving the initial energy 
density p^j in the <fi field, or equivalently the maximum 
temperature of the radiation Tmax- Indeed, one has 
Hi ~ P^j/Mp ~ T^ AX /(Tl n M P ). The latter relation 
can be derived from p<p ~ T 8 /T^ H and the considera- 
tion that the maximum energy in the radiation equals 
the initial (maximum) energy p^j. 

If the neutralino reaches chemical equilibrium, it is 
clear that its final density does not depend on the ini- 



tial conditions. An approximate condition for reach- 
ing chemical equilibrium is [!j (crv) > 10 -9 GeV -2 
(rn x /100GeV)(T RH /MeV)- 2 . Even without reaching 
chemical equilibrium, the neutralino density is insensitive 
to the initial conditions provided the maximum temper- 
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In Fig. 1 we show how the neutralino density fl x h 2 
depends on T RR for illustrative values of the parameter 
r] = 6(lOOTeV/m0), both for neutralinos which are un- 
derdense and which are overdense in usual cosmologies. 

The behavior of the relic density as a function of T RR is 
easy to understand physically. If (av) is large enough so 
that chemical equilibrium is achieved, the usual thermal 
production scenario occurs for T RR > Tf.o.- Neutralino 
annihilation compensates thermal production until the 
latter ceases to be effective at T = Tf. .. The LSP density 
is then determined by the condition r ann = n(av) ~ H 
at T = Tf.o.. Using H ~ T 2 /Mp, as appropriate 
for a radiation dominated Universe, this leads to Yq ~ 
Yto. - (n/a)f.o. - {H/s(av))i. . ~ l/(T f . .Af>(«;)). 
This gives the usual result O s td — {mso/p c )Yo ~ 
(ms )/{p c Tf. o .M P (av)) ~2x 10- 10 GeV" 2 /(fJu). Here 
we used s /p c = 2.8 x 10 8 GeV _1 . When (av) is too 
small for chemical equilibrium to be achieved, the usual 
equation does not hold. Notice that this is the case for 
the smallest (av) in Fig. 1. 

For T R h < 2f. ,, there are four different ways in which 
the density flh 2 depends on T RR : 

(1) Thermal production without chemical equilibrium. 
In this case fi x oc T^ H (e.g. steepest part of the i] = 10" 



(av) = 2 x 10 
estimated in Ref. 
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GeV line). The relic density was 



(av) 



O cdm lO" 16 GeV V 



f lOOGeV y /Trh 
i x J {GeV 



3/2 



(10) 



This matches our numerical calculation for (av) = 2 x 
10- 20 GeV- 2 . 

(2) Thermal production with chemical equilibrium. In 
this case f2 x cx Tp H (e.g. steepest part of the r] = 10~ 7 , 
(av) — 2 x 10~ 12 GeV~ 2 line). The neutralino freezes 
out while the universe is dominated by the <f) field. Its 
freeze-out density is larger than usual, but it is diluted by 
entropy production from (f> decays. The new freeze-out 
temperature is determined by solving n(av) ~ H 
at T = T f N < f w . Using the relations between H, a, and 
T in the decaying-</> dominated Universe, one finds 0, 
13 ^ x - r RH r f.o.(Tf N EW r 4 ^std- Our numerical results 
indicate a slope closer to T^ H , perhaps due to the change 
in Tf.o.. 

(3) Non-thermal production without chemical equilib- 
rium. Here O x cx Trh (e.g. leftmost part of each line). 
Non-thermal production is not compensated by annihi- 
lation. The production of neutralinos is purely non- 
thermal and the relic density depends on rj. It can be 
estimated analytically as follows. For each superpartner 
produced, at least one LSP will remain at the end of a 
chain of decays (due to i?-parity conservation), and thus 
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n x ~ bn^. Here n<f, = p^/ra^. At the time of </>-decay 
p x ~ m x bp^/m^ ~ Tp H , and the entropy is s ~ T^ H . 
Hence p^/s ~ T RH and Y = Y dccay — Wrh/to^- It 
follows that, 

A, 2xl0 3j TO * (ii) 
ficdm V100 GeV/ VMeVy V ; 

(4) Non-thermal production with chemical equilib- 
rium. In this case oc T^" H (e.g. central part of 
the r) = 0.5, (av) = 2 x 1(T 5 GeV" 2 line). Here an- 
nihilation compensates for the non-thermal production 
of neutralinos until the non-thermal production ceases 
at T = T R h- The condition for determining the relic 
density is r ann ~ at T = T RH . This leads to 
Y ~ Y RH ^ IV( Sr hM) ~ l/{T RH M P (<jv)). From 
here it follows that 

fl x ~ (T f . ./T RH )atd. (12) 

With the help of Fig. 1 and formulas (TTJ and {TSJl, 

we can separate the different ranges of 77. Notice that in 
Fig. 1 m x = 100 GeV, thus T f . Q . ~ to x /20 ~ 5 GeV. For 
2rh > Tf. ., the standard production mechanism is re- 
covered, thus Q x = figta (indicated by the horizontal lines 
on the right of the Figure). If T R h < Tf. ., the value of 
f2/i 2 depends on 77. Overdense neutralinos, i.e. those with 
f2 s td > ficdm (above the dotted line labelled WMAP), re- 
quire values of 77 < 10 _4 (100GeV/m x ) to bring their den- 
sity to Cl x — f^cdm- This bound is derived from the BBN 
condition T RR > 5 MeV by taking Q x — f2 c dm m eq. ifTTJl. 
The condition T R h < Tf.o. in Eq. shows that for a 
solution fl x — n c dm with 77 < 10~ 7 (100GeV/m x ) 2 the 
production must be thermal with entropy dilution (case 
(1) or (2)), and T RR must be close to Tf. Q .. Notice that 
in between the two values of 77 just mentioned the pro- 
duction is purely non-thermal (no chemical equilibrium, 



namely case (3)). For all overdense neutralinos, given one 
value of 77 < 10~ 4 (100GeV/m x ) there is only one value 
of T R h for which il x = O c dm- 

Underdense neutralinos, i.e. those with fi st d < ^cdm, 
can have zero, one, or two solutions Q x — f2 c dm- 
There is no solution if Q s td is too low. Neutrali- 
nos with fi s td ^5 10~ 5 (100GeV/m x ) cannot be brought 
to f2 c dm, independently of 77. This can be seen by 
imposing the condition T R h > 5 MeV and fi x = 
J] c d m in Eq. 112(1 . Neutralinos with f2 cdm ^ ^std ^ 
10~ 5 (100GeV/m x ) cannot be brought to f2 cdm ei- 
ther if 77 < 10- 7 (100GeV/77i x ) 2 (fi cdm /f7 std ). This 
happens when the non-thermal production is insuffi- 
cient to increase the density up to r2 c dm, such as 
for the 77 = 10~ 7 , {av) = 2 x 10" 7 GeV~ 2 line 
in Fig. 1. Again using Eqs. (|11[1 an d an d 
Fig. 1, one can see that for the same range of densi- 
ties ricdm ^ ^std ^ 10~ 5 (100GeV/m x ) there are two 
solutions for 10- 7 (100GeV/m x ) 2 (^ cdm /fi std ) < 77 < 
10~ 4 (100GeV/m x ) and a single solution for larger values 
of 77. The two solutions have different values of T RR and 
are both non-thermal one belonging to case (3) and the 
other to case (4). The single solution belongs to case (4). 

In conclusion the neutralino density can be that of cold 
dark matter provided f2 st d <L 10~ 5 (100GeV/m x ) and the 
high energy theory accomodates the combinations of val- 
ues of b/rrirp and T R h identified in the previous two para- 
graphs. These conditions may place constraints on high- 
energy models derived from superstring theories. 
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